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f~^ \ of NRQCD factorization formalism. We find a quite unique feature that the 

0^ \ color-singlet contribution to this process vanishes at not only the leading 

CL|| order but also the next to leading order perturbative QCD calculations and 

Q ' that the dominant contribution comes from the color-octet ^q subprocess. 

The nonperturbative color-octet matrix element oi Sq is related to that of 

kS \ '^S\ of J/ifj by the heavy quark spin symmetry, and the latter can be from 

the direct production of J/ijj at large transverse momentum at the Fermilib 

Tevatron. We then conclude that the measurement of this process may clarify 

the existing conflict between the color-octet prediction and the experimental 

result on the J/^JJ photoprodution. 
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Studies of heavy quarkoniuin production in high energy coUisions may provide impor- 
tant information on both perturbative and nonperturbative QCD. It is well known from 
compelling theoretical |1|J^ and phenomenological [^ reasons that the Color Singlet model 
(CSM) used to describe the production and decay of quarkonium is not complete. In the 
CSM, it is assumed that the cc must be produced in a color-singlet state with the same 
angular-momentum quantum number as the charmonium meson which is eventually ob- 
served. However, it is well know that perturbative QCD calculations of production and 
decay of P-wave quarkonia within the CSM are plagued by infrared divergences 0. More 
recently, measurements made at the Fermilab Tevatron show that the CSM also fails to 
predict the production cross sections of S-wave quarkonia (e.g. J/ip and ip') 0,|^. The 
naive CSM has been supplanted by the NRQCD factorization formalism 0, which allows 
the infrared safe calculation of inclusive charmonium production and decay rates. In this 
approach, the production process is factorized into short and long distance parts, while 
the latter is associated with the nonperturbative matrix elements of four-fermion operators. 
It also predicts a new mechanism called the color-octet mechanism in which a cc pair is 
produced at short distances in a color-octet state, and hadronizes into a final state charmo- 
nium nonperturbatively. This color-octet mechanism might provide an explanation for the 
Tevatron data on the surplus of J/ip and tp' production |PJ^. 

Even though the color-octet mechanism has achieved some successes in describing the 
production and decay of heavy quarkonia, more tests of this mechanism are still needed. 
Recently, the photoproduction data from HERA [§] puts a question about the universality 
of the color-octet matrix elements P,p!0| , in which the fitted values of the matrix elements 



(Cg (^So)) and (Cg (^Pj)) are one order of magnitude smaller than those determined 
from the Tevatron data [|^. (More recently, possible solutions for this problem have been 



suggested in |TT] , p!2| . In this situation it is certainly helpful to find other processes to test 
the color-octet mechanism in the heavy quarkonium production. 

In this paper, we will calculate photoproduction of r^^, 7P — > rjc + X, under the framework 
of NRQCD. The photoproduction processes proceed predominantly through photon-gluon 



fusion at high energies. In NRQCD the Fock state expansion for r]c is 

\Vc) = 0{l)\ccCSo,l)) + Oiv)\ccCPum + Oiv')\cceSu8or l)gg) 

+0{v')\ccCSo,8orl)gg) + ---. (1) 



This expansion has also been given in Ref. ||13[. But, they missed the last term in their ex- 
pansion. This term is crucially important for r]c production in jp process which will be shown 
in the following. For the production of rjc, the contributions to the NRQCD matrix elements 
from the last three terms of the above expansion are the same order of f ^ according to the 
NRQCD velocity scaling rules. They are all suppressed by f ^ compared to the contribution 
from the first term (the color-singlet contribution). However, the color-singlet contributions 
to rjc photoproduction vanish in the leading order and the next to leading order jg fusion 
processes (see the following). So, the photoproduction of rjc is purely a color-octet process, 
even to the next to leading order of QCD calculations. Therefore, the rjc photoproduction 
processes such as at the HERA, will provide an important test for the color-octet produc- 
tion mechanism. Furthermore, we will show by the following calculations, the dominant 
contribution to rjc photoproduction comes from the last term of the Fock state expansion in 
equation (|l|). For this term, the associated color-octet production matrix element {Og^i^So)) 
can be related to the matrix element {Og{^Si)) by the heavy quark symmetry. And the later 
color-octet matrix element is important for the color-octet explanation of the prompt j/ijj 
surplus production at the Tevatron [§|,0. So, the measurement of rj^ photoproduction pro- 
cesses is closely associated with the test of the color-octet gluon fragmentation mechanism 
proposed in p. 

According to the NRQCD factorization formalism [0], the production process 7 + (? — >■ 
T]c + X can be expressed as the following form, 

a{^ + g -^ Vc + X) = J2F{i + g ^ n + X){0:i^). (2) 

n 

Here, n denotes the cc pair configuration in the expansion terms of Eq. (|I]) (including 
angular momentum ^'^"'"^Lj and color index 1 or 8). F{'y + g ^ n + X) is the short distance 



coefficient for the subprocess -^ + g ^ n + X. {O"^) is the long distance non-perturbative 
matrix element which represents the probability of the cc pair in n configuration evolving 
into the physical state rj^. The short distance coefficient F can be calculated by using 
perturbative QCD in expansion of powers of as- The long distance matrix elements are still 
not available from the ffist principles at present. However, the relative importance of the 
contributions from different terms in @) can be estimated by using the NRQCD velocity 
scaling rules. 

From Eq.(|lD, we can see that the color-singlet contribution to the production of rjc is at 
the leading order in v"^. The associated short distance coefficient is given by the subprocess 

^g^cc{^S^,l)+g. (3) 

This process occurs at the next to leading order in a^ for the 75' fusion processes. However, 
there is no contribution from this process, because it violates C(charge) parity conservation. 
(The C parities of the two gluon system (constrained in color- singlet) and the cc pair in 
(^5*0, 1) state are both +1, while the C parity of the photon is —1). 

The color-octet contributions to the rjc photoproduction come from the leading order and 
the next to leading order 75? fusion processes as shown in Fig.l. At the leading order in a^, 
the subprocess is 2 ^ 1 (Fig. 1(a)), 

7^--cc(i5o,8). (4) 

For this process, we readily have 

s 2 
a(7 + <? - cc(iSo, 8) ^ r^,) = -^^5{s - Am^){Of{'So)). (5) 

where s is the invariant mass squared of the partonic process. tjIc is the charm quark mass, 
and we have approximated the charmonium bound state mass of r]c by 2mc. 

At the next to leading order in a^, the subprocesses are 2^2 (Fig. 1(b) and Fig. 1(c)). 
These processes contain the contributions from 



^ + g^ccCSo,8)+g, (6) 

^ + g^ccCSi,8)+g, (7) 

^ + g^ccCPi,8)+g. (8) 

The last two subprocesses only have contributions from diagrams shown in Fig. 1(b). The 
diagrams shown in Fig. 1(c) do not contribute to (^5*1,8) and (^Pi,8) subprocesses because 
of C parity conservation, and they only contribute to the (^5*0,8) production subprocess. 
The cross sections for these 2-^2 subprocesses can be expressed as the following form, 

^(7 + ^ - r^. + X) = j^,FC'^'Lf){Ore'^'Lj)), (9) 

where t = {z — l)s, and z is defined as z = p ■ k^^^/p ■ k^ with p, fc^^, k^ being the momenta 
of the proton, the outgoing rjc and the incident photon respectively. For the short distance 
coefficients F of the subprocesses ^S"} and ^Sq , we readily have P,|l^ 



+ 10sP + 28st^M + 2Qsiu^ + 8su^ + 5t^ + IOPm + 13?m^ + 8tM^ + 4m^). (11) 

Here, s, i and -u satisfy the relation s + i + u = Ami- For the ^Pi subprocess, we calculate 
the short distance coefficient and get 

^n^(8)^ 320(47r)^e^aa^ ,„(.- ^^^ ,„t;^ .^^-^ ,.^.9 „./i-c( 

F(ipf ) = ; ^^ / (g^t + s^M + As'P + Ashu + As'u' + 7sH' 

3[s + ty^[s + uy{t + uY 

+ lls^?M + lls^tM^ + 7s%^ + 7s^t^ + UsH^u + IQsH'^u^ + lAsHv!^ 
+ls^u^ + AsH^ + lls^t^M + IQsH^u^ + IQsH'^u^ + lls^tV + As^iC' 
+sP + Ast'u + llst^M^ + lAsPu^ + IIsPm^ + Asiu^ + sm^ + Pu 
+Ai^u' + Tt^M^ + TFm" + 4t2M^ + tM^) . (12) 

Have provided all the cross section formulas for the subprocesses, we may estimate the 
Tjc production rate in the photoproduction processes. Because there is no color-singlet con- 
tribution in our calculation, all of the following numerical results are for the color-octet 



production. The results depend on the numerical values of a^, rric, and the factorization 
scale Q. We use ^^(m^) = 0.3, nic = 1.486*6^^ and Q^ = (2mc)^. For the parton distribution 
function of the proton, we use the GRV LO parametrization ||15|| . 

In Fig. 2, we show the rjc photoproduction rate via the 2^1 subprocess, in which 
the intermediate color-octet state is ^Sq . Following the heavy quark spin symmetry, we 
estimate the associated color-octet matrix element (O^^i^So)) to be 

{OfCSo)) ^ {OtCSi)) = 1.06 X 10-^GeV\ (13) 

The value of the color-octet matrix element (Og i^^i)) follows the fitted value in [0 by 
comparing the theoretical prediction of direct J/ip production to the experimental data at 
the Tevatron. The 2 — i> 1 process contributes the rjc photoproduction in the forward direction 
(z ~ 1 and p|. ~ 0). In Ref. 0, the authors calculated the J/ip production in the forward 
direction and found some conflicts between the NRQCD prediction and the experimental 
measurements. So, the numerical result shown in Fig. 2 may be also questioned by the same 
problem. However, we note that in the case of J/ip production, the conflicts with experiment 
focus on the values of the color-octet matrix elements {Oq(^So)) and (Og (^Pq))- which are 
determined by the lower pt distribution data of direct J/ip production at the Tevatron ^ . 
The lower px J/'ip production may be strongly affected by the initiate and flnal states gluons 
radiation, which makes the determination of the matrix elements {OglJ'So)) and {0^{^Po)) 
from the Tevatron data very poor |ll[. However, we see that the rjc photoproduction in the 
forward direction of Fig. 2 depends on the value of the color-octet matrix element {O^'' (^ Sq)) , 
which is related to the color-octet matrix element {Og{^Si)) by the heavy quark symmetry. 
And the later color-octet matrix element is determined by the large pt J /ip production data 
at the Tevatron, which can be viewed as a more reliable estimate because in the large pt 
region gluon fragmentation is dominant. 

In Fig. 3-4, we show the photoproduction cross section of rjc via the 2 — > 2 subprocesses 
at the HERA, to which there are three color-octet subprocesses contributions {}Sq , ^S*} 
and ^Pi subprocesses). In Fig. 3, we plot the differential cross section da/dz distribution 
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for Tjc production. In Fig. 4, we plot da/dp'^ distribution. In these two figures, there are 
contributions from ^S*} and ^Pi , for which we estimate the associated color-octet matrix 
elements by using the naive NRQCD velocity scaling rules, 

{OrCS^)) ^ {OfCSo)), (14) 

{orcPi)) ^ (orcso)). (15) 

In these two figures, we impose a cut on z and P^: 0.2 < z < 0.8 for Fig.4 and P^ > IGeV'^ 

1 (8) 

for Fig. 3. From these two figures, we can see that the contribution from the 5*0 subpro- 
cess dominates over the contributions from the other two processes. Furthermore, the ^Sq 
contribution to the z distribution of the differential cross section (Fig. 3) rises rapidly as 
z increases. We note that for the J/ip photoproduction process, the z distribution of the 
cross section from the color-octet contribution has also this property [^. However, this is 
not consistent with the experimental results ||^. This conflict between the color-octet model 
prediction and the experiment invoke a lot of interesting for the further investigations on 
the J/ip photoproduction []TT],0|. In |]TT|, the authors noted that the initiate and flnal state 



gluon radiation may strongly affect the determinations of the color-octet matrix elements 
{Oq(^So)) and {Og{^Pj)) from the direct J/ip production data. This is because the deter- 
minations of these matrix elements are sensitive to the data of J/tp production in the lower 
Pt region, in which the gluon radiation effects are important. These matrix elements are 
also the dominant contributions to the J/ip photoproduction in the high z region. So, by 
involving these effects, they lowered down the determined values of these color-octet matrix 
elements, and found they would lead to a consistent description of the z distribution for 
the J/ip photoproduction. However, these arguments do not affect the rjc photoproduction 
of our results in Fig. 3. This is because the dominant contribution in Fig. 3 comes from the 
color-octet matrix element {O^^I^Sq)), which can be safely related to the matrix element 
{Oq{^Si)) by the heavy quark symmetry, and the value of the matrix element {0^{^Si)) 
we used is from the fltting to the J/ip production data in the large pt [0 region which 



is weakly affected by the initiate and flnal states gluon radiation ||Tl| andtherefore can be 



viewed as a reliable estimate. Another explanation conflict problem emphasizes that the 
intrinsic k^ may be important for the J/ip photoproduction ||12|. If this is true, our predic- 
tions of Fig. 3 and Fig. 4 for the rjc photoproduction will also be affected by the intrinsic kx 
effects. Nevertheless, the measurement of the rjc photoproduction can provide a consistent 
test of these explanations of the J/ip photoproduction conflicts []n],|I2[ . 

For the experimental observation of t]c, one may choose the following decay modes, such 
as T]c -^ 00, rjc — > 77 and rjc — > pp, which either have substantial branching ratios or may 
be easy to detect. From Figs. 2-4, we can see that the production rate of rjc in the photo- 
production processes is comparable with that of J/ip, and we hope that the experimental 
measurement of rjc photoproduction will soon appeare. 

In conclusion, we calculated the rjc photoproduction under the NRQCD factorization 
formalism. We find that the color-singlet contribution to this process vanishes at the leading 
order and the next to leading order QCD calculations. Therefore, this process is a pure color- 
octet process, and the measurement of this process can provide useful information to clarify 
the existing conflict between the color-octet prediction and the experimental results on the 
J/ip photoprodution. 
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Figure Captions 

FIG. 1. The generic Feyninan diagrams for the photoproduction of r]c at the leading order 
(a) and the next to leading order {h)and{c) '-/g fusion processes. All of these diagrams are 
for the color-octet processes. 

FIG. 2. The cross section of 7 + (? ^ ?7c + X in the forward direction {z ~ 1, p^ ^ 0) as 
a function of the photon-proton cm. energy Ec.m. = ^Js-yp- 

FIG. 3. The differential cross sections do jdz for the process 7 -|- p ^ ?7c + AT at the HERA 
as a function of 2; = En^J E^, where ^^s^p = lOOGeV.The carves denote contributions from 
the color-octet ^5*1, ^Pi, and ^5*1 respectively. 

FIG. 4. The differential cross sections da/dP^ for '~f+p —* rj^ + X at the HERA as a function 
oi p^, where ^Js^p = lOOGeV.The carves denote contributions from the color-octet ^Si, ^Pi, 
and ^5*1 respectively. 



11 



7 'WWWi — -~*.^_^^^ 
9 nr6~sir6~ir^ — '*' „;7i 



(a) 



cells'. 



(8)1 



7 'V\AyAA/Vr 



g nrs~5inr5~^ 




ccf^+^if] 



(b) 




(c) 



cc-[i5(«)] 



100- 



T3 
CO 

o 



D 



30 




E_(GeV) 



cm 



100^ 



^o- 



> 

o 

I- 
N 

ID 



1 - 



0.1 




0.2 



0.4 



0.6 



0.8 



(N 



o 

O 
V 

N 

V 



(N 



Oh 



1 - 



0.1 - 




Pt' (GeV') 



